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Star formation in the hosts of high-z QSOs: Evidence from 

Spitzer PAH detections 

D. Lutz 1 , E. Sturm 1 , L.J. Tacconi 1 , E. Valiante 1 , M. Schweitzer 1 H. Netzer 2 , R. Maiolino 3 , 

P. Andreani 4,5 , 0. Shemmer 6 , S. Veilleux 7 



ABSTRACT 



We present Spitzer rest-frame mid-infrared spectroscopy of twelve z~2 mm- 
bright type 1 QSOs, selected from unlensed and lensed QSO samples and covering 
a range of AGN optical luminosities L5100 = 10 45 to 10 47 erg s _1 . On top of 
the AGN continuum, we detect PAH emission from luminous star formation 
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in nine objects individually as well as in the composite spectrum for the full 
sample. PAH luminosity and rest frame far- infrared luminosity correlate and 
extend the similar correlation for lower luminosity local QSOs. This provides 
strong evidence for intense star formation in the hosts of these mm-bright QSOs, 
sometimes exceeding lOOOM yr _1 and dominating their rest frame far-infrared 
emission. The PAH-based limit on star formation rates is lower for luminous 
z~2 QSOs that are not preselected for their mm emission. Partly dependent on 
systematic changes of the AGN dust covering factor and the dust spectral energy 
distribution of the AGN proper, the spectral energy distributions of mm-faint 
high-z QSOs may be AGN dominated out to rest frame far-infrared wavelengths. 
Towards the most luminous high-z QSOs, there is a flattening of the relation 
between star formation and AGN luminosity that is observed for lower redshift 
QSOs. No QSO in our sample has a PAH-measured star formation rate in excess 
of 3000 Moyr- 1 . 

Subject headings: galaxies: active, galaxies: starburst, infrared: galaxies 



Introduction 



Studies of the evolution of galaxies and active galactic nuclei (AGN) are drawing an 
increasingly detailed picture of their global evolution from redshifts above 6 to the present. 
The cosmic star forming density as well as the number density of AGN rise and fall over 
cosmic time, with a maximum in the 'quasar epoch' at redshifts around 2. Less is known 
about the connection of star formation and nuclear activity in individual high redshift ob- 
jects, and the controlling physical processes. A related basic observational goal is to quantify 
both AGN luminosity and star formation in an as wide as possible set of galaxy populations, 
a task that is posing considerable technical difficulties for some types of objects. 

A particularly interesting population are the luminous AGN that are traced by optically 
selected QSOs. Evolutionary connections between intense starburst events and QSOs have 
been proposed early on (e.g. ISanders et al.l Il988l ; iNorman fc Scovilld Il988l ) and naturally 



suggest that the luminous AGN phase coincides with or follows after a period of intense star 
formation. While still limited in spatial resolution and requiring simplifying assumptions, 
increasingly advanced models of these processes now become available. They c ombine galaxy 



merg ing, gas transport, star formation, accretion and feedback processes (e.g. ISpringel et al. 



20051). In turn, such scenarios are part of recen t global modes of galaxy and merger evolution 
(e.g. iGranato et all 120041 : iHopkins et al1l2006h . 
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Indications have been obtained for intense star formation events in the host of high 
redshift QSOs. Perhaps the most important is the detection of submm or mm continuum 
(rest frame far-infrared) from dust in part of the high-z QSO population. If heated by star 
formation, t his emission indicates high star formation rates similar to those in submillimeter 
galaxies (e.g. lOmont et al.ll2003l ; iPriddey et al.ll2003l ; iBarvainis fc Ivisonll2002l ) . Large molec- 
ular gas reser voirs sufficient to sustain intense star formation h ave been detected in CO (see 



summaries m 



Solomon fc Vanden Boutll2005l ; iGreve et al.l 2005 ) , and in some of the brightest 



systems also in tracers of dense molecular gas (IBarvainis et al. 



Carilli et al.l 120051 ; iRiechers et al.l 120061 ; iGracfa-Carpio et al 



1997 



200 



6; 



Solomon et al. 



Guelin et al 



nally, the [CH] 157Atm fine structure line was detected in a few high-z quasars (jMaiolino et al. 



2003 



2007). Fi- 



20051 ; llono et al.l 120061 ) at a ratio to the rest frame far- infrared emission similar to the ratio 
in local ULIRGs, consistent with massive star formation. While these observations indicate 
a plausible picture of massive star formation in some QSO hosts, they need independent 
support, in particular to better discriminate between star formation and AGN heating of 
the observed rest frame far-infrared continuum. We have hence initiated a program using 
the mid-infrared 'PAH' emission features as tracers of star formation in the hosts of high-z 
QSOs. 

Because of their presence in star forming environments with a wide range of physical 
conditions but absence close to powerful AGN, the main mid-infrared PAH emission features 
at 6.2, 7.7, 8.6, and 11.3/im have been used extensively as extragalactic star formation 
indicators. The calibration factor between PAH luminosity and star formation rate depends 
on average physical conditions in the interstellar medium of the galaxy probed. This is 
indicated both in the relatively weaker PAH in galaxies with intense int erstellar radiatio n 
fields (as measured by the large grain dust temperature from their SEDs. iDale et al.ll200ll ). 
and in the sublinear increase of PAH emission with Pa hydrogen re combination line em ission, 
that is observed in spatially resolved studies of nearby galaxies (jCalzetti et al.l 120071 ). For 
the specific case of star formation indicators for QSO hosts, the associated uncertainty is 
outweighed by several advantages. PAH emission from starbursts is very luminous and thus 
less easily outshone by the powerful AGN. Unlike other popular star formation indicators 
like Ha or mid-infrared 24/xm continuum, it will not be quickly overwhelmed by strong AGN 
emission in the very same tracers. It is also less sensitive than e.g. the optical [Oil] 3727A 
line to the significant internal obscuration that is typical for starburst regions. 



In ISchweitzer et al.l ( 120061 ) and iNetzer et al.l (120071 ) we have applied these techniques 
to local z<0.3 Palomar-Green (PG) QSOs, detecting PAH emission in many objects and 
concluding that likely most of their ~ 10 10 -10 12 L Q far-infrared luminosity is due to star 
formation. Comparing average spectra and SEDs of groups of QSOs that are differing in 
their level of far-infrared emission, we were able to derive an infrared SED for the pure AGN, 
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after subtraction of host star formation. 



Several spectroscopic studies have shown that Spitzer-IRS spectroscopy with modest in- 
tegration times of at most a few hours per source can successfully detect the PAH emission of 
the most luminous z~2.5 star forming events, which approach infrared luminosities of 10 13 L Q 
and star formation rates a round lOOOM^/y r . Such detections have been obtained for submil 



limeter galaxies (SMG s) (ILutz et al.l 120051 ; IValiante et al.l 120071 ; iMenendez-Delmestre et al. 



200(1; 


Pope et a 


2008 


) as well as for 


Weedman et al. 


2006|; 


Yan et al. 


2007) 



in an IRS spectral order have been demonstrated to be able to detect yet fainter PAH 
emission at similar redshifts with a peak of the rest frame 7.7/zm feature of 0.4mJy or less 
( ITeplitz et al.l 120071 ). Probing PAH emission to these levels would be highly interesting for 
a complete sample of luminous z~2.5 QSOs but expensive in observing time. Systematic 
difficulties arise in detecting PAH features on a strong AGN continuum that additionally has 
structure due to silicate emission. The silicate emission is strong in regions of the spectrum 
where the main PAH complexes are weak, and thus tends t o wash out the PAH-in duced 
structure, as demonstrated for some local PG QSO spectra in ISchweitzer et al.l (120061 ) . 



In this paper, we hence focus as an important first step on measuring PAH emission 
and inferred star formation in a modest size sample of twelve z=1.8 to z=2.8 QSOs that 
all have detected (sub)mm emission. If this mm emission is due to host star formation, its 
associated PAH emission will be det ectable with IRS in modest integration times by analogy 
to the SMGs at similar redshift. In lLutz et al.l (120071 ) we have presented first results for the 
brightest object in our sample, the Cloverleaf (H1413+117). Here we present the full mm- 
bright sample and discuss in section 4.2 implications for the high redshift QSO 1 population 

whole. 

Throughout the paper we adopt an Q m = 0.3, = 0.7 and H = 70 kms _1 Mpc _1 
cosmology. When converting IR lum inosities to star fo rmation rates we adopt for consistency 
with much of the previous work the iKennicuttl (119981 ) value which is based on a O.l-IOOM© 
Salpeter IMF. For a more realistic Chabrier or Kroupa IMF star formation rates would be 
a factor ~1.6 smaller. Our sample includes lensed sources. For these, we present fluxes or 
flux densities as observed, but luminosities, star formation rates and masses corrected for 
the adopted lensing magnification. 
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2. Sample and Observations 



We have selected from the literature z~2 Type 1 radio quiet QSOs with robust (> 5a) 
and consequently bright submm or mm dust continuum detections. The sample includes 
lensed as well as unlensed QSOs and is summarized i n Table [p Sou r ces ha ve be en chosen from 



the su rveys of (sub)mm emission in z^ 2 QSOs by Omont et al.l ( 20031 ) and IPriddey et al. 



Guilloteau et al. 



( 2003), the lensed QSO sub mm survey of Baryainis fc Ivison ( 20021 ). and one target each from 



(119991) and lKnudsen et al.l (120031 ) . One of the targets (SMMJ04315+10277, 



Knudsen et al.l (120031 )) was first reported in a submm survey and may be intermediate in 



properties between optically selected QSOs and typical submm selected galaxies. The other 
targets are selected from the mentioned studies of optically selected QSOs solely on the basis 
of a bright and robust (sub)mm detection and will provide a view of the rest frame far-infrared 
bright part of the general z~2 unobscured QSO population. Our unlensed targets have AGN 
rest frame optical luminosities L5100 = vL u (5100A) in the range 10 46 to 10 47 ergs s _1 , some 
of the lensed QSOs are intrinsically an order of magnitude less lu minous and overlap with 
the highest luminosity l ocal sources studied with similar methods by ISchweitzer et al.l (120061 ) 
and iNetzer et al.l toO^ \. 



For targets in lOmont et al.l (120031 ) the optical luminosities were adopted from the rest 
frame blue absolute magnitudes in this reference, converti ng to our adopted co smology. In 
the other cases we used published near-infrared photometry ( IPriddey et al.ll2003l . and 2MASS 
point source catalog), and the assumption of a f u oc z/ -0 44 rest frame optical/UV continuum 
JVanden Berk et alJlSoOlh . 



We obtained low resolution ( R~ 60 — 120) mid- infrared spectra of the QSOs using the 
Spitzer infrared spectrograph IRS (IHouck et al.ll2004j ) during the period July 2006 - October 
2007. The IRS LL slit is wide enough to include all images of the lensed sources. For each 
target, 30 cycles of 120 sec integration time per nod position were taken in the LL1 (19.5 to 
38.0 fim) module and 15 cycles in the LL2 (14.0 to 21.3 /xm) module, leading to effective on- 
source integration times of 2 and 1 hours, respectively. We start from the pipeline processed 
basic calibrated data (version 14.4 to 16.1 for different sources), subtracting individually 
data for the same data collection event within the data collection sequences of the two nod 
positions. We use own deglitching routines identifying outliers from the default bad pixel 
masks, by identifying points discrepant by more than a multiple of the local noise from their 
spectral neighborhood, and by visual inspection of the 2D frames, finally replacing them with 
values representative of their spectral neighborhood. We coadded the images for each cycle, 
checking pixel per pixel for outlying values in the stack, and subtrac ted residual backgr ound 
derived from a first order spatial fit to source-free regions. SMART ( iHigdon et al.ll2004l ) and 
SPICE (optimum extraction, for the weak continuum source SMMJ04315 only) were used 
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for extraction. We compared the independent spectra from the two nod positions before 
averaging, finding no mismatches that would be indicating mispointings or other serious 
problems with the data. Before combining the segments of the spectrum observed in the 
different modules into the final spectrum, the segment observed in the LL2 module was 
multiplied with an object-dependent scaling factor (within a few % of 1) in order to provide 
best match in the overlapping wavelength range. 



3. Results 



Figure [TJ shows the spectra of all our targets, and Fig. [3] a composite derived from all 
sample objects. The main spectral features are the ones also seen in the spectra of local 



PG QSOs fe.g. ISiebenmorgen et al.ll2005l ; lHao et al.ll2005l ; ISchweitzer et al.ll2006l ; lHao et al. 



20071 ; iNetzer et al.l 120071 ): A strong AGN continuum with superposed silicate emission at 
A > 9/iin, and aromatic 'PAH' emission features at 6.2, 7.7, 8.6/irn at varying strength 
relative to the local continuum. An unusual individual spectrum is the one of SBS 1408+567 
with convex continuum, a smooth broad emission feature around 8/xm, but no unambiguous 
identification of PAH emission. This broad feature is weaker relative to the continuum but 
has similarity with the 8/xm fe ature in the mid-infrared spectrum of the well known Type 



1 AGN and ULIRG Mr k 231 flRigopoulou et al.l 1199a ISpoon et al.l 12002 : IWeedman et al. 



20051 : lArmus et al.l 120071 ) which may originate from a combination of silicate emission and 
absorption. PAH 7.7/im may also contribute to some extent to this feature. 

In this paper, we focus on analysis of the PAH emission and its implications for star 
formation in the QSO hosts. Given the S/N of the spectra and the variations in slope and 
shapes of continuum and silicate emission in AGN, we do not attempt to derive independent 
fluxes for the various PAH features. Informed by the spectra of local QSOs and supported by 
the average spectrum of Fig. [31 we quantify the overall level of PAH emission by decomposing 
the QSO spectra into a scaled and resampled starburst template (based on the ISO-SWS 
spectrum of M82. ISturm et al.ll2000l ) at the nominal redshift of the QSO, and the underlying 
AGN emission which consists of a dust continuum and silicate emission that is turning up 
beyond ~8/im. Results of this decomposition are shown in Figure [2j The combination 
of AGN continuum and silicate emission varies smoothly with wavelength over the range 
of interest. The first signature of PAH in such spectra that is already visible at modest 
S/N then is the strongest 7.7/im feature and its short wavelength steep flank, while its 
longer wavelength side is often less pronounced, due to the combined presence of 8.6/im 
PAH and the upturn of silicate emission. The strength of the starburst component in the 
decomposition is adjusted to produce a smooth residual AGN emission at the locations of 
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the 6.2 and 7.7/^m features, in particular to remove the steep short wavelength flank of the 
7.7/im feature. Given the S/N of the data we do not attempt to derive PAH feature ratios by 
multi-feature fits, but obtain from the simple template a more robust characterisation of the 
overall level of PAH emission which we represent by the flux of the brightest (7.7/xm) feature. 
These PAH 7.7/xm fluxes (Table [2]) are measured in the scaled starburst template resulting 
from the decomposition, by fitting Lorentzian features for the PAHs plus an underlying 
residual continuum. PAH 7.7/xm fluxes listed in Table [2 ] can t hus be directly compared to 



the results for local QSOs presented by ISchweitzer et al.l (120061 ). 



We list in Table[2]the measured noise per resolution element over the rest frame 6.5-8/xm 
region of the actual spectra, defined as RMS around a linear fit to the data after subtraction 
of the scaled PAH template and its associated feature structure. By quoting the noise 
measured in this way from the spectra, we implicitly include possible noise contributions for 
stronger continua, that may arise for example due to more difficult outlier identification and 
fringing. This noise may be conservatively overestimated if the continuum is slightly curved 
over this range. Noise ~0.2mJy per spectral sample as obtained for some of our weaker 
continuum targets is in good agreement with the nominal IRS sensitivity for our observing 
parameters, as implemented in the Spitzer Science Center SPEC-PET tool. 

For the 9 sources for which we report PAH detections in Table [2J the peak height of 
PAH emission in the starburst component obtained from the decomposition is required to be 
at least 3 times the measured local noise. These significances of 3-13 for the PAH emission 
effectively reflect the significance of the short wavelength flank of the 7.7fim feature. At 
longer wavelengths, systematic variations in shape of underlying continuum and in particu- 
lar the strength and detailed shape of the onset of AGN silicate emission at > 8/xm start to 
become more and more limiting, preventing reliable use of these wavelengths to verify or im- 
prove the formal significance. Several of the targets (SMMJ04315+102177, KUV0806+4037, 
HS1002+4400, LBQS1230+1627B) show PAH with peak height that is for the individual 
sources just above the peak/noise=3 detection limit. In combination with the possibility 
that the mentioned uncertainties on AGN continuum and silicate shape affect to a lesser ex- 
tent also the short flank of the 7.7/im PAH, we estimate that individual PAH fluxes for these 
weak PAH detections can be uncertain at the 30-50% level. We caution that they should 
not be used for analyses going significantly beyond the modest accuracy requirements of the 
luminosity analyses presented below. 

Upper limits are estimated analogously, for the PAH template scaled to a level where the 
features would remain visibly detectable on top of the continuum and the silicate feature. The 
IRS spectra of the sources with PAH nondetections are still consistent with the presence of 
PAH emission at peak feature fluxes similar to the sources with detections. Again, the main 
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limitation is systematic, due to the strong and structured underlying continuum and silicate 
emission in these targets. In the spectrum of LBQS 0018-0220, for example, inspection of 
Fig. [T] shows that the 7-8/xm region of the spectrum, observed with good S/N, may represent 
either PAH emission partly filling in the region before an upturn of silicate emission, or be 
due to variations in shape/slope of this feature, which is strong for LBQS 0018-0220, and 



the underlying continuum. Such variation s are observed in local AGN (e.g. lHao et al.l 12005 



Schweitzer et al.ll2006l ; iNetzer et al.ll2007l ). We list in Table [2] PAH upper limits considering 
this possiblity rather than only reflecting the detector noise. 

We will compare below the PAH emission of high-z QSOs with the rest frame far-infrared 
emission. Partly, this is possible by simply using the rest frame ~ 250/zm continuum that 
can be obtained from the published (sub) mm photometry with minimal extrapolation. For 
other purposes and for comparison with the local QSO population we need to extrapolate 
to the rest frame far-infrared peak which is constrained by direct observations for very few 
high redshift QSOs only. We do so by assuming that the detected (sub)mm emission reflects 
a rest frame T=40K, emissivity index (3=1.5 greybody spectrum. These parameters are 
in reasonable agreeme nt with the dust temp erature ~36K suggested for the s ubmillimeter 



galaxy population by I Chapman et al.l (120051 ) and with 40-60K suggested by iBeelen et al. 



(120061 ) for high-z QSOs, from the combination of 350/mi photometry with longer wave- 
length observations. For thr ee QSOs ove r lappin g with our sample (KUV 0806+4037, HS 
1002+4400, SBS 1408+567). IBeelen et al.l kood ) find 32K, 38K, and 35K respectively, and 
the FIR-warmer QSOs in their small sample may include objects with a smaller ratio of host 
FIR/submm to AGN emission. Similarly, our simple extrapolation to rest frame 60jj,m for 
H14 13+117 (the Cloyerleaf) is in good agreement with the observations of iBenfordl (119991 ) 
and iRowan- Robinson! (120001 ) . At redshift z~2.4, changing the adopted temperature by +5K 
will change the estimated rest frame 60/im flux by +60% and —45%, respectively. This is 
the dominant source of uncertainty in the estimated rest frame 60/im emission, since the 
original (sub)mm measurements for our IRS sample are S/N>5 by sample definition. Using 
the T=40K (3=1.5 assumption we extrapolate to the rest frame FIR peak of this component 
fo r all our sources, which we express as uL u (Q0fj,m) for consistency with the local QSO work 
of ISchweitzer et al.l (120061 ). These 60/im luminosities are listed in Table [2J For sources with 
published independent 850/zm and 1.2mm photometry we have averaged results from the 
two possible extrapolations. 
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4. Discussion 

4.1. Consistent evidence for luminous starbursts from PAH and rest frame 

far-infrared 



Comparison of the PAH detections and limits to the rest frame far-infrared emission of 
the sample QSOs, and comparison to other populations is the key to testing the hypothesis 
that both PAH and far-infrared trace intense starbursts in the hosts. We do such comparisons 
first in directly observed quantities, and then using PAH luminosities in combination with 
FIR luminosities derived with the greybody extrapolation described above. 

Fig.H]shows a comparison of the PAH 7.7/im peak flux density with the flux at rest frame 
wavelength 222/im. The latter corresponds to an observed wavelength of 850Mm for z~2.8 



the specific redshift chosen for consistency with previous work on SMGs (ILutz et al.l 12005 



Valiante et al.ll2007l ). This is a simple diagram based on directly observed quantities, and has 
minimal sensitivity to SED assumptions for the extrapolation from the observed submm/mm 
point to rest frame 222/zm. The distribution of the ratio of continuum-subtracted PAH 
peak flux density and rest frame FIR continuum Sp Am. 7/ W ' S^ jim is very simil ar for our 



mm-bright QSOs, for the submillimeter galaxy sample of IValiante et al.l (120071 ). and for 
local ULIRGs. The median log(SpAH7.7^m/ is -1.14, -1.32 and -1.26 for the QSOs, 
SMGs, and ULIRGs, respectively, and the distribution for the QSOs can be drawn at >25% 
probability from the same distribution as either the SMGs or the ULIRGs, according to a 
two-sided Kolomogorov-Smirnov test. 

Moving from observed quantities to PAH and far- infrared luminosities, Fig. [5] shows that 
the mm-brig ht high-z QSOs e x tend the relation between these luminosities found for local 
PG QSOs by ISchweitzer et al. J2006h by an order of magnitude towards higher luminosities, 
and that the relation for the QSO population crosses the locus for star-forming ULIRGs. The 
PAH and FIR data of the mm bright QSOs agree with a scenario where these two tracers are 
jointly originating in starbursts in the hosts. These starbursts reach SMG-like luminosities 
of ~ 1O 13 L in the most luminous (unlensed) targets of our sample. 

Individual PAH-based star formation rates can be derived for our sample and are listed 
in Table [21 Concerning the calibration factor applied, the similarity of the PAH/FIR ratio 
in our sample objects and ULIRGs per mits us to use the rati o L(PAH)/L(FIR)~ 0.013 



derived for the starburst-like ULIRGs in Schweitzer et al 



(2006), in combination with the 



calibration of FIR-based star formation rates by iKennicuttl (Il998l ). Because of the changes 
of PAH/FIR with galaxy properties discussed above, it is worth recalling that the resulting 
SFR [Mq year' 1 } = 3.5 x 1Q~ 42 Lpah [ er 9 assumes ULIRG-like physical conditions 
in the interstellar medium of the galaxy observed. PAH based star formation rates for 
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our high-z sample range from about 200 to 3000 M ^/yr. As was found for local QSOs 
and their far-infrared emission (jSchweitzer et al.ll2006l ). use of the PAH diagnostic provides 
strong support for massive star forming events in the mm bright QSOs, at yet an order of 
magnitude larger luminosities. 

Molecular gas CO emission has previously been searched for in six of our targets, de- 
tecting five (Table El all ga s masses based on the conversion factor 0.8 M /(K km s" 1 pc 2 ) of 
Downes fc Solomon! (119981) and our adopted c osmology). Another target (HE 1104-1805) has 



been searched for CO (IBarvainis et al. 



2002]), but the nondetection may not be meaningful 



given the narrow bandwidth of this CO search, which may miss the target CO line for plau- 
sible uncertainties of the optical input redshift. Consistent with their intense star formation, 
the QSO hosts are gas-rich systems (molecular gas mass of a few 1O 1O M ). Despite the large 
gas masses, gas exhaustion times M(H 2 )/SFR are short and range from 10 7 to 8 x 10 7 years, 
where measured. 



4.2. Prom mm-bright QSOls to the unbiased QSOl population 



The previous section has discussed properties of the mm-bright part of the z~2 luminous 
QSOl population only. As an indication of the significance of this s election, we note tha t we 
actually observed 6 m m-bright QSOs out of the combined samples of lOmont et all (120031 ) and 
Priddey et all (120031 ) . which in total contain 83 distinct objects. We discuss the properties of 
o ur objects in the context to these two samples, the 1.7<z<3 objects among the lensed sample 
of IBarvainis fc Ivisonl ( 120021 ) . and the local PG QSOs, in order to elucidate implications for 
the complete and unbiased population of luminous z~2 QSO Is. 

Figure [6] shows the QSO PAH luminosity as a function of AGN rest frame optical 
luminosity. The mm bright Q SOs clearly extend t o higher luminosity the relation for local 
PG QSOs already discussed by lNetzer et all (120071 ). with a possible saturation (flattening of 
the relation) at the largest AGN luminosities. An extension to individual PAH detections or 
significant limits for the full population including mm-faint QSOs would be extremely difficult 
even with Spitzer capabilities. Substantial insight can already be gaine d, nevertheless, by 
adding to Fig. [6] the limit on PAH emission set by iMaiolino et all (120071 ) for the average of 
a similar redshift QSOl sample that is not mm prese lected. Unlike their average spectrum, 
the individual spectra from the IMaiolino et all (120071 ) sample are of insufficient S/N to put 
a limit on PAH emission that is meaningful for our analysis. Compared to mm bright QSOs 
from our sample of similar AGN luminosities, the limit of PAH emission for the not mm- 
biased average spectrum is clearly lower, by a factor ~3 compared to our PAH detections. 
The corresponding average SFR limit is <7OOM yr~ x . 
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In an analogous way, Fig. [7| compares the QSO far-infrared luminosity, for which a larger 
number of individual measurements is available, to the AGN rest frame optical luminosity. 
This diagram also reminds that the z~2 mm-bright QSOs, having L FIR ~L 510 o i-e. L FIR ~0.1 
\jboU are far from being extreme infrared excess objects and resemble in this property the 
average local QSO. In relation to the AGN luminosity, the strongest IR emission for the 
high-z QSOs in this diagram is for SMMJ04315+1027, which was originally detected at 
submm wavelengths. In this diagram, the bulk of the SMG population would extend from 
the location of SMMJ04315+1027 to the left, because of the typically modest lumin osity 



of the AGN they are often hosting (e.g. Alexander et al.l 12005 



Valiante et all 120071 ). As 



concerns the relation of mm bright QSOs to the total QSO population, we obtain a similar 
result as from the PAH data. Again, the mm-bright QSO studied here extend the relation 
for local QSOs and represent the far-infrared bright end of a distribution in L^/^/L 510 o 
among the z~2 QSOs. Considering only the L 510 o > 10 46 erg s _1 QSOs, the mm bright QSOs 
have median LpvR=10 46 ' 6 ergs _1 , a factor of 4 larger than the median for the full population. 
The specific value of this factor should be considered tentative only, since the distribution 
for the full population is based on (sub)mm photometric points of typically low individual 
significance. 



m 



The PAH spectroscop y pre sented for mm brig ht high-z QSOs here, and for local QSOs 



Schweitzer et al.l (120061 ) and iNetzer et al.l (120071 ). has provided strong support for a star 



formation origin of the rest frame far-infrared emission in these systems. This conclusion has 
been derived for QSOs with L FIR ~L 510 o, however, and needs to be critically reviewed for 
th e even more AGN d ominated objects at the mm-faint end of the high-z QSO population. 
In INetzer et al.l (120071 ). we have presented a host-subtracted intrinsic AGN SED for local 
PG QSOs which is very broadly characterised by z/L y (5100A)~ z/L^/im)^ 5 x z/L,,(60/ihi), 
with considerable scatter between individual objects. This intrinsic SED drops significantly 
towards the rest frame far-infrared, but in combination with the properties of the high-z 
QSOs (Fig. [7]) would imply that the mm-faint end of the z~2 QSO population in fact may 
contain almost 'pure' AGN without strong host star formation. The AGN-heated dust could 
already provide the observed weak rest frame far-infrared emission. A caveat to be noted here 
originates from recent evidence for systematic changes in the mid-infrared t o optical SEDs 



of luminous high-z QSO s compared to their local low- luminosity analogs (jMaiolino et al. 



2007 ; iTreister et al.ll2008l ). and from the possible detection of high-z AGN without hot dust 
(I Jiang et al.l 120061 ) . We defer a detailed analysis of the AGN-heated dust in our sample 
to a future paper, but note that the larger PAH equivalent width of our average spectrum 
compared to the local PG average (Fig. [3]) is in support of these results, given that it is 
combined with similar or lower ratio of PAH and 5100A continuum in our sample compared 
to locally (Fig. [6]). These results suggest a lower mid- infrared AGN continuum luminosity 
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in relation to the optical AGN luminosity, and thus smaller AGN dust covering factors. 
They imply c hange s of the intrinsic AGN SED for luminous high-z QSOs compared to the 



Netzer et al.l (120071 ) intrinsic SED for local QSOs, likely weakening the AGN-related 60/im 



emission along with the weakening of the rest frame mid-infrared AGN dust emission. This 
decrease of AGN rest frame 60/im emission would reduce again the part of the total far- 
infrared continuum directly ascribed to the AGN and perhaps suggest a noticeable host 
contribution even in mm-faint high-z QSOs. Firmly establishing the nature of the weak 
rest frame far-infrared emission at the mm-faint end of the z~2 QSO population will hence 
require a more complete SED characterisation over the full rest frame infrared range. 



4.3. A flattening relation of starburst and AGN luminosity 



For local QSOs with L5100 ~ 10 44 to 10 46 ergs _1 , we have found clear correlations of 
the star formation indicators PAH and far-infrared continuum with AGN optical luminosity 
( jNetzer et al.l 120071 ). that are approximately linear over this luminosity range. The high-z 
QSOs extend the AGN luminosity range by another order of magnitude. The relation of 
PAH (star formation) and AGN luminosity is flattening at L5100 > 10 46 ergs _1 (Fig. [6]). This 
is not yet unambiguous from the comparison of our mm-bright sample alone with the local 
PG QSOs. There may be a slight flattening, and the slope of a linear fit in Log(Lp ah) vs. 
log(L 5 ioo) is 1.01 for the PAH detections among the local PG QSOs, versus 0.92 adding the 
high-z PAH detections to the sample, but these changes are not yet significant. The effect 
becomes clear by incl usion of the PAH upper limit which has been derived for the sample of 
Maiolino et al.l (120071 ) that is not selected to be mm-bright (Fig. [6]). In agreement with mm 
continuum based studies, our PAH spectroscopy does not find evidence for star formation 
rates above 3OOOM yr _1 (starburst luminosities much larger than 1O 13 L ). This saturation 
occurs at luminosities similar to the upper end of the luminosity range spanned by the submm 
galaxy population, a limit that is often related to the star formation rate in a maximal 
starburst that is co nverting much of its gas into stars over just a few dynamical times (e.g. 
Tacconi et al.l 120061 ) . In support of this view, a slower than linear increase (or saturation at 
high AGN luminosities) is also inferred for the relation of cold g as mass and AGN lumi nosity, 
from a synopsis of CO measurements for high-z/high-L QSOs (jMaiolino et al.ll2007bl ). 



Local PG QSOs are found in a region of the fundamental plane for elliptical galaxies and 
spheroids that is hosting moderate mass ellipticals with corre sponding black hole masses ~ 
10 8 M Q , and accreting efficiently at an Eddington ratio ~0.25 (IDasyra et al.ll2007l ). Relative 
to the host star formation rate, their AGN growth rates are more than an order of magnitude 
larger than needed to explain the locally observed BH mass to bulge mass relation under 
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the assumption of coeval growth (INetzer et al.ll2007l ). The increased L(AGN)/L(SF) for the 
high-z QSOs and in particular the mm-faint ones places them in a yet more pronounced 
phase of AGN growth proceeding much more rapidly than host growth. The failure up to 
now to identify matching extreme >> 1O 14 L star forming events also in other galaxies and 
surveys indicates this is not just due to a time delay between AGN and host growth, but 
must reflect intrinsically different durations of the AGN growth and host growth phases. 



5. Conclusions 

In the context of quantitatively constraining the coexistence of star formation and AGN 
activity in different classes of high redshift galaxies, we have used Spitzer-IRS mid-infrared 
spectroscopy to detect PAH emission in the spectra of nine out of twelve mm-bright z~2 
radio quiet QSOs, with optical luminosities L 5 i o = 10 45 to 10 47 erg s -1 . The detections 
provide strong support for the presence of intense star formation up to ^OOOMQyr" 1 in 
their hosts, that is dominating their rest frame far-infrared emission. The typical luminous 
z~2 QSO that is not preselected for mm emission has lower star formation rate <7OOM yr _1 
and the mm-faint part of this population may include objects that are AGN dominated out 
to rest-frame far-infrared wavelengths, due to luminous AGN but at best modest level of 
star formation. Toward the high luminosity and redshift QSOs, there is a flattening in the 
relation between starburst and AGN luminosity that is observed for local PG QSOs, and no 
object in our sample has a PAH-measured star formation rate above 3OOOM yr _1 . 

This work is based on observations made with the Spitzer Space Telescope, which is 
operated by the Jet Propulsion Laboratory, California Institute of Technology, under a con- 
tract with NASA. Support for this work was provided by NASA under contracts 1287653 
and 1287740 (S.V.,O.S.). 
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Table 1. Sample 



Name 


z 


RA 


DEC 


fx L Ref. log(L 5 ioo) 






J2000.0 


J2000.0 


erg/s 



LBQS 0018-0220 


2 


596 





21 


27 


30 


-2 


03 


33 





1 







46 


82 


HE 0230-2130 


2 


162 


2 


32 


33 


10 


-21 


17 


26 





14 


5 


B02 


45 


29 


SMMJ04315+10277 


2 


837 


4 


13 


27 


26 


+10 


27 


40 


5 


1 


3 


K03 


45 


70 


KUV 0806+4037 


1 


795 


8 


12 


00 


54 


+40 


28 


13 





1 







46 


30 


RXJ0911.4+0551 


2 


793 


9 


11 


27 


61 


+5 


50 


54 


1 


20 





B02 


45 


32 


HS 1002+4400 


2 


097 


10 


05 


17 


43 


+43 


46 


09 


3 


1 







46 


84 


HE 1104-1805 


2 


305 


11 


06 


33 


40 


-18 


21 


23 





10 


8 


B02 


45 


70 


LBQS 1230+1627B 


2 


735 


12 


33 


10 


40 


+ 16 


10 


52 





1 







46 


89 


SBS 1408+567 


2 


583 


14 


09 


55 


57 


+56 


28 


26 


5 


1 







46 


89 


H1413+117 


2 


558 


14 


15 


46 


27 


+ 11 


29 


43 


4 


11 





V03 


46 


11 


HS1611+4719 


2 


35 


16 


12 


39 


90 


+47 


11 


57 





1 







46 


61 


J164914.9+530316 


2 


26 


16 


49 


14 


90 


+53 


03 


16 





1 







46 


89 



References. 
K03: 



Lensing magnifi cations i ij, are adopted from: B02: iBarvainis fc Ivison 



Knudsen et al.l (120031 ) . V03: IVenturini fc Solomon! (12003( 1 . Throughout the 



((20021) 

paper, we present fluxes or flux densities as observed, but luminosities, masses, and star 
formation rates corrected for these adopted magnifications. 
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Table 2. Results and host starburst properties 



Name 


Peak (PAH7.7) 
mjy 


Noise a 
mjy 


F (PAH7.7) 
10~ 21 Wcm~ 2 


L (PAH7.7) 
10 44 erg s — 1 


fZ/„ (60/m?i) 
10 44 erg s" 1 


SFR (PAH) 

Meyr- 1 


M (H 2 ) 
1O 1O M 


LBQS 0018-0220 


<1.9 


0.27 


<2.3 


<12.0 


543 


<4200 


<3.5 b 


HE 0230-2130 


1.7 


0.21 


2.3 


0.6 


48 


197 




SMMJ04315+10277 


0.8 


0.20 


0.8 


4.8 


615 


1680 


13. 3 b 


KUV 0806+4037 


1.1 


0.30 


1.7 


3.7 


392 


1300 




RXJ0911.4+0551 


2.7 


0.21 


3.0 


1.0 


43 


340 


0.5 b 


HS 1002+4400 


1.1 


0.35 


1.6 


5.0 


337 


1760 




HE 1104-1805 


2.8 


0.36 


3.8 


l.i 


43 


496 




LBQS 1230+1627B 


1.2 


0.28 


1.1 


8.1 


609 


2840 


2.4 C 


SBS 1408+567 


<3.0 


0.32 


<3.6 


<19.6 


888 


<6840 


6.6 b 


H1413+117 


5.1 


0.75 


6.1 


2.9 


117 


1020 


4.0 d 


HS1611+4719 


1.0 


0.19 


1.2 


5.3 


394 


1840 




J164914.9+530316 


<1.5 


0.24 


<2.0 


< 7.9 


399 


<2780 





a Measured noise per spectral sample in the rest frame 6.5-8^tm range 
lHainline etal] ll2004ll 
iGuilloteau et alj <199Sft 
< 1Weifi et al.l 11200311 
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Fig. 1. — Spitzer-IRS spectra of the sample of (sub)mm-bright high-z QSOls. Vertical lines 
indicate the location of the 6.2 and 7.7/im PAH features at the redshift of the sources with 
PAH detections. 
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Fig. 2. — Decomposition of IRS spectra. The individual panels are shown in the same 
ordering and scaling as in Fig. [H with the thin black lines retracing the spectra as observed. 
Thick blue lines show the AGN continuum and silicate emission remaining after subtraction 
of the PAH component obtained in the decomposition. The green line indicates the la noise 
measured in IRS data obtained with same observing set up and integra t ion ti me but lacking 
a source (one of the nondetections among the SMGs of IValiante et al.l (120071 )). While noise 
from this curve in good agreement with expected IRS sensitivity and measured noise for 
our continuum-weaker QSOs and gives an indication of trends with wavelength, it misses 
possible sources of noise in spectra with bright continuum that can arise e.g. due to more 
difficult identification of outlying values. Instead, we conservatively present in Table [2] noise 
measured in the PAH region of the spectra after subtraction of the PAH template. This 
noise is the RMS after subtracting a linear fit to the rest wavelength 6.5-8/im region (before 
the main upturn of silicate emission). 
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Fig. 3. — Average rest frame mid-IR spectrum of the sample high-z QSOls. The individual 
spectra have been scaled before averaging to the same flux of 15mJy at rest frame wavelength 
222/mi. This allows a direct comparison of spectrum and the relation of PAH and fa r-infrared 
emission to the stacked submm galaxy spectrum presented by lValiante et al.l (120071 ). A direct 
average of the QSO s pectra has a very sim ilar shape. A scaled average spectrum of the local 
PG QSOs studied by iNetzer et al.l (120071 ) is shown for comparison. 
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Fig. 4. — Ratio of PAH peak flux density and rest frame 222/im far-infrared continuum for 
the sample QSOs (top panel). z~2 SMGs and loc al ULIRGs are added for comparison in 
the lower panels (SMGs and ULIRGs adopted from Ivaliante et alll2007h . 
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Fig. 5. — Relation of the PAH 7.7/xm luminosity and far-infrared continuum luminosities for 
our sample and for local comparison objects. P G QSOs as well a s the compa rison sample 



of sta rbur sting ULIRGs are from the studies of ISchweitzer et ajj (120061 ) and iNetzer et al. 



(120071 ) . For lensed high-z objects, luminosities have been corrected for the magnifications 
listed in Table [H 
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Fig. 6. — Top panel: Relation of PAH 7.7/im luminosity and AGN 5100A continuum lumi- 
nosity. The large 'A verage high-z QSO' upper limit symbol refers to the limit derived by 



Maiolino et al.l (120071 ) for the average spectrum of a high redshift QSO sample which is not 



selected for bright (sub)mm emission. Bottom panel: To help identifying trends, the same 
data are shown plotting the ratio of PAH 7.7/im luminosity and AGN 5100A continuum 
luminosity as a function of AGN 5100A continuum luminosity. 
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Fig. 7. — Top panel: Relation of 60/xm continuum luminosity and AGN 5100A continuum lu- 
minosity. The 60/iin luminosities of high redshift QSOs are extrapolated from their (sub)mm 
photometry using a T=40K (3=1.5 greybody approximation. The green hexagon indicates 
the locus for the median of the full z~2 L 5 i o > lO^ergs" 1 population. Bottom panel: 
To help identifying trends, the same data are shown plotting the ratio of FIR luminosity 
and AGN 5 100 A continuum luminosity as a function of AGN 5 100 A c ontinuum luminosity. 
' Othe r' hig h-z QSOs refers to o bjects in the studies of lOmont et al.l (120031 ). iPriddey et al. 
(120031 ). and iBarvainis fc Ivisonl (120021 ) that are not part of our mm-bright sample. 



